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ABSTRACT 


The maximum fuel pin temperature for nominal operation in an advanced 
power reactor is 1370 K. Because of possible nitrogen embrittlement of the 
clad, the fuel temperature was limited to 1622 K. Assuming simultaneous 
occurrence of the most adverse conditions a "deterministic" analysis gave a 
maximum fuel temperature of 1610 K. A statistical analysis, using a syn- 
thesized estimate of the standard deviation for the highest fuel pin tem- 
perature, showed probabilities of 0.015 of that pin exceeding the tempera- 
ture limit by the distribution free Chebyshev inequality and virtually nil 
assuming a normal distribution. The latter assumption gives a 1463 K max- 
imum temperature at 3 standard deviations, the usually assumed cutoff. 
Further, the distribution and standard deviation of the fuel-clad gap are 
the most significant contributions to the uncertainty in the fuel tempera- 
ture. 
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PRELIMINARY ANALYSIS OF HOT SPOT FACTORS IN AN ADVANCED 
REACTOR FOR SPACE ELECTRIC POWER SYSTEMS 
by Paul H. Lustig, Arthur G. Holms, and Harry W. Davison 
Lewis Research Center 


SUMMARY 

Hot spot factors have been calculated for the fuel pins in a lithium 
cooled Advanced Space Power Reactor. These hot spot factors were estim- 
ated by calculating maximum temperatures in the uranium mononitride fuel. 
They account for uncertainties in power, coolant flow, material proper- 
ties, and fuel pin dimensions. The hot spot factor is the ratio of the 
difference between the maximum fuel temperature and the coolant inlet 
temperature, with the uncertainties included, to the temperature differ- 
ence for the nominal design condition. 

The maximum fuel pin temperature for nominal design conditions is 
1370 K (2006° F) . A fuel temperature limit of 1622 K (2460° F) was set 
to prevent possible nitrogen embrittlement of the clad. This is a pre- 
liminary and, supposedly, conservative limit. 

A "deterministic" analysis indicated a maximum temperature only 12 K 
(21.6° F) below an assumed fuel temperature limit of 1622 K (2460° F) . 

This marginal safety factor under the assumption of simultaneous occur- 
rences of the most adverse conditions indicates that a statistical analy- 
sis is desirable to obtain a more representative assessment of reliability. 

A less pessimistic estimate of the maximum fuel pin temperature is 
1463 K (2173° F) arrived at through a statistical analysis which assumes 
a normal distribution of the fuel temperature and a cutoff at three stand- 
ard deviations. The probability of exceeding this temperature is, then, 
0.0013. Using another statistical approach, an estimate was made of the 
probability of exceeding the fuel temperature limit based on the Chebyshev 
inequality which does not assume any particular mathematical function for 
the frequency distribution of the fuel temperature. This probability 
could be as much as 0.015 for the central pins. 

The results of this study indicate the probability distribution func- 
tion of the fuel-clad gap to be the most significant contributor to 
the uncertainty in fuel temperature. 


INTRODUCTION 

The performance of a reactor fuel pin is determined in part by the tem 
peratures occurring within it. It is desirable to know the distribution of 
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temperature not only for design conditions , but also for any deviations 
from design which may possibly appear in an actual fabricated fuel pin 
functioning in a reactor core. Such information can be used to help guide 
the design by dictating materials, shapes, and dimensions, and setting re- 
quirements for tolerances. For instance the analysis can be used. to es- 
tablish which tolerances must be close and which may be relaxed from the 
point of view of fuel temperature effects. 

This report deals with off-design effects on the maximum fuel pin 
temperature in an Advanced Power Reactor (APR) for Space Electric Power 
Systems investigated at the Lewis Research Center. It is desirable to 
know how much the temperature of a given spot may be influenced by off- 
design conditions. In particular, it is desirable to know what the maxi- 
mum temperature in the fuel pin may be. Deviations from their design 
values may occur in a number of factors influencing temperature, and sev- 
eral such deviations may occur in the same pin, at the same place in the 
pin, and at the same time. However, it is not likely that all factors 
will have maximum deviations at the same location and at the same time. 

Two approaches are used to estimate the maximum fuel pin temperature 
in the off-design condition. One approach, which is very conservative, 
assumes that all the variables are at their extreme expected values coin- 
cidentally and that each operates on the fuel pin temperature to produce 
its maximum value. This is the deterministic method which in its general 
form amounts to a worst case analysis. Amendola (ref. 1) formulated a 
deterministic method as a product of individual hot channel factors ap- 
proximated by an abbreviated Taylor series expansion. 

The. second, more realistic, approach assumes the variables to be 
random statistical quantities. It is possible then to determine the fre- 
quency distribution of the fuel pin temperature by a method described by 
Abernathy (ref. 2). An alternative to that employs the assumption 

of a linear combination of independent variables but this time to obtain 
an approximation of the standard deviation of the fuel pin temperature. 

The standard deviation is needed to calculate the probability of exceed- 
ing a fuel temperature limit and to calculate the hot spot factor. The 
standard deviations of the independent variables are based on tolerances 
and published data since there is not yet any data on sample variations 
for this particular design. 

In this report, the hot spot factor is calculated using both deter- 
ministic and statistical approaches. The hot spot factor is the ratio 
of the temperature difference between maximum fuel temperature and inlet 
coolant temperature for the off -design conditions and the temperature 
difference for the reference conditions, respectively. One technique 
(ref. 3) for quoting the hot spot factor by a statistical method is to 
assume that the fuel temperature for the off-design conditions is a tem- 
perature such that is has a probability of only 0.13 percent of being 
exceeded. This corresponds to three standard deviations of the fuel 
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temperature. However, as will be discussed later, the preliminary nature 
of the design dictates that a conservative estimate of probability of ex- 
ceeding a fuel temperature limit be made and for this the Chebyshev in- 
equality will be employed. Both the normal distribution approximation and 
Chebyshev inequality are probability estimation methods but the former as- 
sumes a particular distribution of the fuel pin temperature while the 
latter makes no assumption of a, particular form of the frequency distri- 
bution. 

This report includes a description of the fuel pin design, a discus- 
sion of the various uncertainties in variables affecting temperature, and 
a summary of the analytical techniques employed. The assumptions of in- 
dependence and linearity of the variables are discussed briefly. The de- 
sign radial and axial temperature distributions are discussed in Appendix B 
A typical calculation of the hot spot factor is given in Appendix C. 


DESCRIPTION OF REACTOR AND FUEL PIN 

The reactor being studied has a fast spectrum and is fueled with 
uranium mononitride. An isometric drawing of the reactor is shown in fig- 
ure 1. Each of 181 fuel pins is placed within a T-lll honeycomb structure 
to maintain a triangular pitch. The honeycomb structure is surrounded by 
an annular neutron reflector of TZM (Mo-0.5Ti-0.08Zr). The entire struc- 
ture is surrounded by a T-lll pressure vessel which is 58 » 5 centimeters 
(23 in.) in diameter and 68.6 centimeters (27 in.) long. Reactivity is 
controlled by rotating six TZM drums in the annular reflector. Eleven 
fuel pins are placed in holes drilled on one side of each drum and an an- 
nular segment of the drum opposite the fuel pins contains a T-lll neutron 
absorber. 

A sketch of a typical fuel pin and coolant channel are shown in fig- 
ure 2. The UN fuel is clad with 0.147 centimeter (0.058 in.) thick T-lll. 
A 0.0127 centimeter (0.005 in.) thick tungsten liner physically separates 
the fuel from the T-lll clad. This liner inhibits any chemical reaction 
between the UN fuel and tantalum in the clad. A clearance provided be- 
tween the fuel and liner for ease of assembly, is filled with helium. Dur 
ing initial operation this nominal radial gap width is about 0o0038 cm, 
and closes as the fuel swells. Fission product gas collection spaces are 
provided at both ends of the pin and at the center of the fuel (the fuel 
cylinder is hollow). Fuel spacers are located at both ends of the fuel 
column to allow fuel swelling in the axial direction and to inhibit vibra- 
tion of the fuel during launch conditions. Support buttons, located 120 
degrees apart and at five axial positions on the honeycomb support tube, 
center the fuel element within the coolant channel and inhibit bowing. 

The reactor design power is 2.17MW and the axial power shape can be 
approximated by a chopped cosine curve having a peak to average power, 
ratio of 1.23. The radial power distribution across the reactor is flat- 
tened to 1.33 by using 3 fuel zones. The desired fuel loading in each of 
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the three zones is obtained by adjusting the diameter of the central hole 
in the fuel (table I). Pins containing the most fuel are located fur- 
thest from the center of the reactor. The greatest power density occurs 
in the center pin at the reactor midplane. 

The fuel is cooled by liquid lithium which enters the reactor at 
1165 K (1640° F) . The primary lithium flow is in an annular space be- 
tween the fuel and the honecomb . The average coolant temperature rise 
across the reactor is 55 K (100° F) . The convective film coefficient is 
large so the clad wall temperature is nearly equal to the coolant temper- 
ature. The core design radial and axial temperature distributions are 
discussed in Appendix B. 


DISCUSSION OF UNCERTAINTIES 

Ten variables which contribute to the value of the highest fuel tem- 
perature in the reactor have been considered. The variation of each in 
the following description is quoted at two standard deviations except for 
the variables that are strictly a result of tolerance limits. These ex-, 
ceptions are local fuel thickness, orifice mismatch, flow annulus width, 
fuel-clad gap, fuel shift in the clad, and clad shift in the flow annulus. 
For these, the number of standard deviations associated with the half- 
range H is /3 or 1.732. 

The entire list of variables is as follows: 

(1) Overall Power Adjustment. Reactor power may be determined from 
measurements of fluid flow rate and coolant temperature difference from 
inlet to outlet. Variations from the desired reactor power can result 
from inaccuracy in the measurement of flow and temperature and time lags 
in the response pattern of the controls. 

An uncertainty of +8 percent on overall power adjustment was assumed 
(see ref„: 3). Although this variation (+8 percent) is associated 

with three standard deviations in reference 9, in this study the ±8 per- 
cent uncertainty is associated with only two standard deviations for con- 
servatism. 

(2) Local Power Density. The power density at a point in the fuel depends 
on the fuel density, enrichment, relative neutron flux, and neutron en- 
ergy spectrum at that point.. The fuel density and enrichment depend on 
the fuel fabrication methods while the relative neutron flux and spectrum 
depend on the nuclear properties and geometry of the materials in the re- 
actor.. An uncertainty of ±4 percent in local power density is assumed. 

(3) Local Fuel Thickness. The estimated uncertainty in fuel thick- 
ness is ±0.0038 centimeter (±0.0015 in.). This uncertainty is based on 
fuel fabrication tolerances. 
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(4) Local Thermal Conductivity of the Fuel. The thermal conductivity 
of the fuel depends on the fuel porosity and temperature. The uncertainty 
in thermal conductivity, estimated from the data of references 4 and 5 , is 
±20 percent. 

(5) Orifice Mismatch. Uncertainty in coolant flow results from ina- 
bility to size orifices to compensate accurately for fabrication discrep- 
ancies in flow passages and in the distribution headers. Tolerances sug- 
gest a value of ±5 percent. 

(6) Flow Annulus Width. Here again the uncertainty is based on fab- 
rication tolerances. These allow a variation of. +0.0121 centimeter 
(±0.00475 in.) about the design value of 0.10 centimeter (0.040 in.). 

(7) Heat Transfer Coefficient. Reference 6 shows data for liquid 
metal heat transfer at Peclet numbers below critical with Nusselt numbers 
varying from 4.9 to 7.8. This is roughly a range of +25 percent about an 
average. On this basis, we assume an uncertainty of +25 percent. 

(8) Fuel-Clad Gap Width. Tolerances in fuel outer radius, tungsten 
liner thickness , and clad inner radius give gap widths at temperature 
varying from 0.0013 centimeter (0.0005 in.) to 0.0063 centimeter (0.0025 
in.). The uncertainty is therefore +0.0025 centimeter (+0.0010 in.). 

(9) Fuel Shift in the Clad. This could occur as a result of fuel 
rod bowing, and cause the fuel-clad gap width to be non-uniform. The 
amount of non-uniformity can be stated in terms of local increase or de- 
crease in width. The value of ±0.0038 centimeter (±0.0015 in.) was chosen 
which is the maximum gap within tolerances for the hot fuel assembly. The 
uncertainty amounts to +100 percent of the gap and includes all possibili- 
ties except, perhaps, such effects as out of round distortion. 

(10) Clad Shifts in Flow Annulus. This depends on fabrication toler- 
ance and also on warpage due to differential expansion resulting from cir- 
cumferential temperature variations in fuel pin and annulus wall. The 
distortion is limited by the dimples in the wall. Uncertainty is based on 
the variation within this limited range. In terms of maximum local in- 
crease or decrease in. annulus width, this amounts to ±0.020 centimeter 
(+0.008 in.). 


SUMMARY OF ANALYTICAL TECHNIQUE 

The maximum fuel pin temperature in the reactor occurs at the center- 
line of the center pin slightly above core midplane. Its nominal value 
Tf m is a function of all the variables that affect it directly when they 
have their nominal values. 

* ***** * 

^fm = T ( x i ,x 2 ,x 3 ,x 4 ,X 5 ’ ’ ' " ,x n^ 


( 1 ) 
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The maximum fuel temperature can be expressed as follows: 


T fra * T<x l + a . x l- x 2 + Ax 2' x 3 + Ax 3 * 


• > x n + 4x a> (2) 


If all the increments are taken at such extreme values as to pro- 
duce maximum temperature , the result corresponds to the general treat- 
ment of the deterministic approach or worst case analysis. As an ana- 
lytical simplification, the temperature deviation or increment Tf m - Tf m 
can be approximated as the sum of the linear terms of the Taylor series 
expansion (ref. 2). 
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3x, 


3T,. 3T r 

. , fm A . fm . . 

Ax, H — Ax„ + — Ax, + 
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— Ax 

3x n 
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(3) 


Two important assumptions are implied in the use of the above equa- 
tion. One is that the variables X 2 ,X 2 » . . . , x n are mathematically 
independent, that is, there are no functional relations among the vari- 
ables. In addition, the relation between the maximum fuel pin tempera- 
ture and the , other variables must be linear or nearly so in the range 
in which the equation is applied. 

The older technique (ref „ 1) for calculating Tf m - T* m had been 
the simple multiplication of the individual hot spot , factors with the 
Axi's at their maximum possible value. This technique is called the 
deterministic method. These factors, described in reference 1, can be 
expressed as follows : 


change in T^ due to uncertainty "i" 
T* - T 

fm coolant inlet 


(4) 


The numerator of the fraction in equation, (4) can be calculated from 
equation (3) if the variable has a nearly linear effect on temperature 
and if the coefficients 3Tj m /3x£ are known. Otherwise, the temperature 
change can be calculated from an expression relating the maximum fuel pin 
temperature to a variable with every other variable at its nominal value . 
In either event, a conservative result is likely to be obtained when the 

n 

product of individual hot spot factors, namely, I I f^ is calculated for 

i=l 

the overall hot spot factor. In fact, if the input variables are inde- 
pendent and the temperature is linear with respect to the input variables, 
then the "deterministic" method yields the same maximum temperature as a 
worst case analysis. 


in contrast with the deterministic approach , a technique has recently 
evolved (ref. 7) that applies some statistical theory of random variables 
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to the input uncertainties and their effect on the maximum fuel pin tem- 
perature. If equation (3) is used with its attendant assumptions then 
the standard deviation of the maximum fuel pin temperature 

ax can be expressed in terms of the standard deviations, a^'s of the 
input variables. Thus 



4 + 





(ref. 2) 


(5) 


Hence, if the standard deviations of the input variables can be estimated, 
then it is possible to estimate the standard deviation of the maximum fuel 
pin temperature. No assumptions have been made at this point as to the 
frequency distribution of the input variables (xf's). 

Once the standard deviation of the maximum fuel pin temperature has 
been calculated, it is possible to estimate the probability of exceeding 
some temperature limit only by making some assumptions about the frequency 
distribution of the fuel pin temperature. 

At this point, most of the literature makes one or the other of the 
following assumptions (in addition to those underlying equation (5)): 

(1) the input variables each have a normal frequency distribution; hence 
the output variable Tf m also has a normal distribution, or (2) the 
Central Limit Theorem applies as an approximation to the author's case. 

The Central Limit Theorem states that if equation (3) is applicable then 
the distribution of Tf m will approach the normal distribution as n 
approaches infinity independently of the distributions of the x^'s. The 
validity of the second assumption depends on the number of variables that 
contribute to the major portion of the variance of Tf m . Thus, Abernathy 
(ref. 2) shows that with only five independent variables whose variances 

contribute 97y percent of the output variance and are within a factor of 

8 of each other, the output variable "closely approximates" a normal dis- 
tribution, even though the input variables had rectangular distributions. 
However, as will be discussed later, the approximation is not necessarily 
close for small or "tail" probabilities. 

Finding the distribution of the output variable may also be accom- 
plished directly by applying the method described by Abernathy (ref. 2). 
This would use the "guessed at" frequency distributions of the input var- 
iables. In general, the method would be tedious and the result, at this 
preliminary stage of design, of such questionable application, that such 
an approach would not be worth the effort. 

Another alternative which yields a very conservative limit as the 
probability function of the output variable based solely on the standard 
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deviation calculated from equation (5) is the application of the Chebyshev 
inequality (ref. 7). Use of Chebyshev 's inequality requires only that 
the variable have a mean and standard deviation. The inequality can be 
stated in the following way as applied here: The probability that the 

maximum fuel pin temperature will exceed Tf m + noq> is less than 1/n^ 
where n is the number of standard deviations, Tf m is the population 
mean of the distribution, and o^ is the value of the standard deviation. 
Thus, if T^ is the fuel temperature limit then the probability of ex- 
ceeding T^ is - 1 /((Tl - T"f m )/0T)2- Assuming that Tf m and ox are 
known exactly, Chebyshev 's inequality is the most conservative approach 
that could be taken at this point since it assumes nothing with respect 
to the frequency function of the temperature distribution except the ex- 
istence of a mean and variance. In view of the preliminary nature of the 
design, it was considered that an estimate of reliability should be made 
based on the Chebyshev inequality. In addition, the probability, based 
on the assumption that the frequency distribution of the maximum fuel pin 
temperature is a normal one, will be calculated. This is the usual pro- 
cedure from the literature and is certainly much more optimistic than the 
Chebyshev inequality probability. It has some validity in the fact that 
the Central Limit Theorem does provide that a linear combination of a 
large number of independent random variables does tend toward a normal 
distribution. Once the frequency distributions are more . precisely defined 
then it would be possible (in theory) to completely define the distribu- 
tion of the maximum fuel pin temperature by a technique such as the one 
described by Abernathy (ref. 2). 


THE ASSUMPTIONS OF INDEPENDENCE AND LINEARITY 

Essential to the development of equation (5) are the assumptions of 
independence among the input variables and linearity between the output 
and the input variables . 

Of the ten input variables considered in this analysis (see table II) , 
the only question as to any dependency arises with respect to those vari- 
ables which affect the coolant flow rate, and the heat transfer coefficient. 
Variables x^ (orifice mismatch) , xg (flow annulus width) , and xjq (clad 
shift in flow annulus) all affect coolant velocity which in turn can alter 
the coolant heat transfer coefficient. However, the critical Peclet number 
for lithium is 300 (ref. 6) and the Peclet number for the coolant. in the 
flow annulus is about 70. This means that the heat is transferred from the 
fuel to the coolant by molecular conduction and is essentially independent 
of the fluid velocity. Hence, the heat transfer coefficient is independent 
of those variables affecting fluid velocity. However, it was found that 
the clad shift in the flow annulus does cause a change in the average heat 
transfer coefficient making the latter a dependent variable. But the in- 
fluence, of the uncertainties in the heat transfer coefficient on the maxi- 
mum fuel pin . temperature is so small that the effect of any interdependence 
on the standard deviation estimate or the worst case analysis can be 
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neglected. There is no dependence that is in evidence among any of the 
other factors. 

The linearity of the relation between the fuel temperature and the 
input variables is assumed without proof. A more detailed investigation 
of the degree of non-linearity and the error introduced by any non- 
linearity is left for a future investigation. Some discussions of Taylor 
expansions to be used in the presence of dependencies and non-linearities 
are given in references 7 and 80 


CALCULAT IONAL PROCEDURE 

The calculation is made for the reactor centerline fuel pin except 
where otherwise stated. Heat generation and temperature in this pin are 
circumferentially uniform if fuel thickness and density are uniform, if the 
fuel remains centered in the clad, and if the clad remains centered in 
the coolant passage. For this ideal situation, since we choose 'to make 
the hot spot analysis at the plane where the fuel temperature is a maxi- 
mum, a radial distribution in the pin is all that is necessary and refer- 
ence conditions are found with a one-dimensional calculation. 

However, the fuel is not, in general, apt to be centered in the clad, 
nor the clad in the coolant passage. Hence, a two-dimensional calculation 
is required. This has been done using the steady state heat transfer pro- 
gram STHTP. The program uses a rectangular nodal scheme for handling one, 
two, or three dimensional problems with heat generation, constant or tem- 
perature dependent thermal conductivity, constant contact or film coeffi- 
cients, and radiative heat transfer with constant emissivity. Cylindrical 
nodes can be handled with sufficient accuracy if radius ratios are not too 
large. 

Node . arrangement for the calculations are shown in figure 3 for the 
reference case where fuel is centered in the clad and clad is centered in 
the coolant passage. Temperature drop across the liner was included in 
the calculation but was insignificant. 

Figure 4 shows the arrangement for the calculations with asymmetry in the 
gas gap and in the coolant passage. 

The hot spot calculation is confined to finding a hot spot factor at 
the beginning of core life assuming 100 percent helium in the gap between 
the fuel and the clad. Even though fission products escaping into the 
helium gap. during irradiation decrease the thermal conductivity of the 
gas during operation, the growth. of the fuel decreases the gap. As a re- 
sult, the maximum fuel temperature at the beginning of core life will be 
greater than any other time during operation. 

The procedure for finding the effect on temperature of the various 
independent factors is, first to make temperature calculations using a 
reference set of conditions. This establishes T of equation (1) where 
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the reference conditions are the nominal values of the variables. Then 
values of expected maximum variation for each factor affecting tempera- 
ture are established. The method and values were discussed earlier in 
the section Discussion of Uncertainties. These are the Ax's of equa- 
tions (2) and (3). New temperatures are then calculated using the max- 
imum perturbed value of the factor or variable. A AT, the difference 
in temperature between the reference case and perturbed case, is then 
calculated. 

Even though we have established what were termed maximum variations 
for each factor, we must establish a degree of confidence that the ac- 
tual value will fall within those maxima. This degree of confidence is 
related to the standard deviation for a particular distribution. To each 
assumed maximum variation, a particular number of standard deviations is 
assigned. The less confidence there was in the maximum variation, the 
fewer the number of standard deviations were assigned to it, which therefore 
produced relatively larger estimates of the standard deviations. 

The particular number of standard deviations assigned to the expected 
maximum variation of each variable was discussed in the section on uncer- 
tainties. For those variables that were considered to have approximately 
rectangular distributions, the number of standard deviations is H/a = /T 
where H is the half range of the variable (ref. 1) . These include the 
variables concerned with dimensional tolerances. Estimation of the stand- 
ard deviation of the remaining variables are derived primarily from the 
literature. There was an element of uncertainity as to the applicability 
of the literature estimates to our own case. Where this element was 
strong, the degree of confidence in the value was reduced by decreasing 
the number of standard deviations quoted. 

The temperature of the fuel pin is calculated for each factor at its 
perturbed value. Each difference between nominal and perturbed tempera- 
tures (Tf m £ -Tf m ) is then divided by the number of standard deviations es- 
timated for that factor. This means all the temperature differences now 
correspond . to one standard deviation in each factor. This AT^ is then 
the term (BT^/Sx^Ax^ in equation (3), with the Ax{s equal to one stand- 
ard deviation. Also when squared, it becomes the term (9Tf /9xi) 2 cr| in 
equation (5) so that the sum of all such squared terms is the 

variance of the maximum fuel pin temperature. This means equation (5) can 
be rewritten as follows when a linear relation exists between input and 
output variable : 

2 2 2 2 

o T = AT, + AT, + » • -AT (6) 

113 n 

where the AT's are the temperature rises accompanying a one standard de- 
viation perturbation of the input variables. 

The hot spot factors and probabilities are to be calculated based on 
assumptions concerning the frequency distribution of the maximum fuel pin 
temperature. The usual technique is to assume that the distribution is 
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a normal one and that at three standard deviations, the probability of 
exceeding that temperature (0.0013) is sufficiently small as to call it, 
highly unlikely. It is with this three standard deviation temperature 
that the hot spot factor is calculated. 

As described earlier, the "deterministic" method was used as a worst 
case analysis approximation. The hot spot factor resulting from this 
technique is also quoted. The hot spot calculations using both the sta- 
tistical method and deterministic method are given in Appendix C. 

Apart from the hot spot factors, a calculation was made for the prob- 
ability of exceeding the maximum fuel temperature limit. of 1622 K (2460 F) . 
This limit is based on a preliminary estimate established to prevent pos- 
sible nitrogen embrittlement of the T-lll clad. This limit was established 
in connection with the work described in reference 9. 


CALCULATIONAL INPUT 

Data required for making the calculation included information on node 
connections and dimensions, as stated in the preceeding section. Also in- 
cluded are thermal conductivities, either constant or functions of temper- 
ature, emissivities for thermal radiating surfaces, convective heat trans- 
fer coefficients, heat generation rate, and sink temperatures. Since the 
calculation is made on a pin cross-section at one core location, the cool- 
ant temperature at that location is taken as the sink temperature. Volu- 
metric heating rate is assumed to be independent of fuel radius and no 
heat is generated elsewhere. Convective heat transfer is based on Dwyer's 
correlation (ref 0 10)= 

Lithium properties are those reported by Davison (ref. 11). Proper- 
ties of the T-lll cladding were obtained from reference 12. UN properties 
are from references 4 and 5 and helium properties from reference 13. 


RESULTS 

The values of the perturbed temperature differences that correspond 
to one standard deviation of the contributing factor or variable are sum- 
marized in table II » The components of the variance or (ATi)^ terms of 
equation (6) are also listed. The magnitude of one standard deviation of 
the maximum fuel pin temperature is 31.0 K (55.8° F) . Applying the re- 
cent conventions of using three standard deviations in the calculation of 
the hot spot factor the result is a fuel hot spot factor of 1.46. The 
corresponding clad hot spot factor is 1.59 as shown in table III. 

The "deterministic" method of multiplying individual hot spot fac- 
tors was also used. The extreme unfavorable value of each variable had 
to be selected. Since a "deterministic" analysis should be a fully 
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pessimistic approach as far as the coincidence of unfortunate events is 
concerned, the additional pessimism of assuming each variable's maximum 
value to be at three standard deviations was used. This was done even 
though the value on a tolerance limited variable was well outside its 
physical possibility. This, however, could compensate for gross errors 
in inspection or some relaxation of tolerances as the design and fabri- 
cation proceeds. The individual hot spot factors are listed in table V. 
The cumulative hot spot factor is 2.18 so that the maximum fuel pin tem- 
perature under "deterministic" (approximate worst case) conditions would 
be 1610 K (2438° F) . 

In the design study, a conservative fuel temperature limit of 1622 K 
(2460° F) has been established in order to avoid the possibility of em- 
brittlement of the clad due to the formation of nitrides from the nitro- 
gen vapor in the fuel. This means that the "deterministic" analysis in- 
dicates only a small safety margin of 12 K (21.6° F) between the maximum 
fuel pin temperature and the fuel temperature limit. This situation war- 
rants a closer look at a more realistic approach inherent in a statisti- 
cal hot spot factor analysis. One such method gives a maximum fuel pin 
temperature of 1463 K using the assumptions that the frequency distribu- 
tion is normal and that the probability of exceeding three standard de- 
viations is essentially zero (for exact normality and exactly known mean 
and variance it is 0.13 percent). 

The reactor will thus be safe with respect to the fuel temperature 
limit of 1622 K by this analysis. This limit stands 8.1 standard devia- 
tions away from the nominal fuel pin temperature. However, as illustrated 
on page 235 of reference 7, while the normal distribution is often a good 
approximation to a real situation in the central range of a distribution 
(say out to two standard deviations), it becomes increasingly in error at 
the tails (namely beyond two standard deviations) . The assumption of nor- 
mality can be avoided using the Chebyshev (or Tchebychev) inequality. ^ 
Using it, the probability is 0.015 that the central fuel pin will exceed 
the fuel temperature limit of 1622 K. 

This represents a small probability that the central fuel pin will 
exceed the fuel temperature limit. However, there are 247 pins in the 
reactor, albeit, at lower temperatures. Let Pp^ be the probability of 
failure or the probability that the fuel temperature limit will be ex- 
ceeded by the i^ pin of the N or 247 pins. The success probability 
for each pin P s ^ = (1 - Ppi) was calculated based on its maximum opera- 
ting temperature. 


l-In using the Chebyshev inequality, it is necessary to assume that 
the nominal or design temperature Tf m coincides with the mean of the . 
temperature distribution T fm- This is a reasonable assumption in view 
of the pessimism implied in the use of the Chebyshev inequality and that 
highly skewed distributions are not expected for any of the variables. 
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The probability that all the pins will operate below the limit is 
N 

the product ~[T P si* assumes the events are independent of each 

1 

other from pin to pin. Most, of the variables are local with respect to 
each fuel pin. Only two of the variables, namely, overall power and heat 
transfer coefficient are global with respect to all the fuel pins. Fur- 
thermore, the fuel clad gap width greatly overshadows all the other vari- 
ables in importance in contributing to fuel pin temperature variability 
(see table III). In view of these considerations, the 247 different fuel 
pins will be assumed to survive or fail, statistically independently of 
each other. The probability that at least one pin fail is simply 
N 

1 ~ TT ^si“ 

1 


The maximum nominal operating temperature was calculated for each 
pin and the standard deviation scaled to that for the hottest pin in pro- 
portion to the amount of heat generated. With these numbers, the Chebyshev 
inequality was used to calculate the failure probability for each 

pin. The probability of failure of at least one pin at the core could be 
as much as 0.82. In other words, until the standard deviations of the in- 
put variables can be better defined and reduced and/or the frequency dis- 
tribution of these variables or the output variable (fuel temperature) are 
known or can be assumed with more confidence, fuel pin integrity with re- 
spect to the failure of the clad by nitrogen embrittlement . should remain 
somewhat in doubt. It should be emphasized that pessimism upon pessimism 
has been heaped on this calculation as far as expanded standard devia- 
tions, worst possible frequency distributions, and the fuel pin tempera- 
ture limit itself is concerned. As indicated previously, even the "de- 
terministic" analysis only shows a marginal chance of the central pin fail- 
ing by this mode, hence, less chance by any other pin. 

Another point that should be made is that while the overall power ad- 
justment is the second greatest contributor to the uncertainty in the fuel 
pin temperature, it is only a transient effect and would probably not have 
any great influence on failure by nitrogen embrittlement of the clad which 
is a long term effect. However, elimination of this uncertainty reduces 
the standard deviation by only 1 K. 

The calculations showed that the fuel-clad gap is the most signifi- 
cant factor in determining the fuel hot spot factor. In this respect, the 
dimension should be kept strictly within tolerance. The uncertainty of 
this factor represents 28.2 K (50.8° F) of the total calculated standard 
deviation of 31.0 K (55.8° F)(see table II). The frequency distribution 
of the fuel-clad gap will be the prime influence on the distribution of 
the fuel temperature. This means that since the frequency distribution 
of the gap can probably be described as somewhat rectangular in shape, 
use of the Chebyshev inequality, where no assumption is made as to the 
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distribution, or in a sense assumes the worst possible distribution, is 
too pessimistic even at this point. 

The variables, local fuel thickness, orifice mismatch, heat trans- 
fer coefficient, and flow annulus space have little effect on the maxi- 
mum fuel temperature. Also, any shifting of the fuel in the clad or 
shifting of the clad in the flow channel would cause a much smaller per- 
turbation in maximum fuel temperature (less than 2 K (3° F) for one 
standard deviation) . 


CONCLUDING REMARKS 

Calculations have been made to find temperature distributions in 
some of the hottest fuel pins in a reactor core, and to estimate the prob- 
ability of extreme temperatures to be encountered in the fuel and in the 
clad. Hot spot results are preliminary because of the preliminary nature 
of the design information used, the paucity of material property data, and 
the lack of hardware for use in estimating condition and behavior of the 
parts involved. The following results should, however, be useful in the 
continuing design work: 

1. The nominal maximum fuel temperature occurs in the centerline 
fuel pin and has a value of 1370 K (2006° F) . 

2. Temperature decreases radially, the maximum nominal values in 
Zones II for fuel being 1353 K (1982° F) and in Zone III is 1333 K 
(1950° F). 

3. A "deterministic" analysis of the fuel pin temperature yields a 
maximum temperature of 1610 K (2438° F) . This indicates only a small mar- 
gin of safety below the 1622 K (2460° F) which it is estimated would cause 
failure due to nitrogen embrittlement of the clad. 

4. A recently favored statistical approach gives a maximum fuel pin 
temperature of 1463 K (2173° F) . The technique assumes a normal frequency 
distribution of the fuel temperature and sets three standard deviations 
above the mean or nominal value as the upper fuel temperature that can be 
expected. The standard deviation in fuel temperature by this technique is 
31.0 K (55.8° F) . The fuel pin hot spot factor is 1.46. 

5. The probability of exceeding a fuel temperature limit of 1622 K 
(2460° F) in the central fuel pin could be as much as 0.015. This state- 
ment is based on the Chebyshev inequality and assumes only that the mean 
and standard deviation of the fuel temperature are known. Thus, there is 
no assumption as to the distribution as was made in the calculation of 
the maximum expected fuel temperature in 4 above. The method is extremely 
pessimistic, but does show the need for a reasonable grasp on the fre- 
quency distributions of the variables. For if we extend the method to all 
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247 fuel pins, the probability for the failure of one or more of them by 
exceeding the limit could be as much as 0.82 which is much too high and 
could occur if the fuel temperture has a very unfavorable frequency dis- 
tribution. 

6. In this light, the probability distribution function of the 
fuel-clad gap is the most significant contributor to the 

standard deviation of the fuel pin temperature. Hence, the statis- 
tical nature of this variable should be closely studied to insure the in- 
tegrity of the fuel pins. Even at this point, it is possible to say that 
the frequency distribution of this variable will be somewhat rectangular 
in shape making probabilities based on the Chebyshev inequality too 
pessimistic. 
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APPENDIX A 
SYMBOLS 

cross sectional area 

arbitrary constants 

fuel temperature rise 
hot spot factor 

half range of symmetrically distributed variable 

thermal conductivity 

number of standard deviations 

probability of failure 

probability of success 

heat generation rate per unit volume 

heat transfer per unit area 

radius 

temperature 

fuel temperature limit 
temperature difference 
ratio of fuel volume to cell volume 
fuel-clad gap width 

variable which contributes to fuel temperature uncertainties 
change in variable from its mean value 

temperature rise from location j to location j + 1 when 
perturbed by off-design condition of i^ variable 

temperature rise from location j to location j + 1 

standard deviation 

standard deviation of dependent variable 
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Subscripts : 
cl clad 

f fuel 

i variable (i = 1, 2, 3, . . 10) 

in inner 

j location along path of temperature rise (j = 1, 2, . . .,6) 

m value of variable at the position where it achieves its maximum value 

til 

n n or last variable 

ref reference case 
Superscripts: 

* nominal or design value 


denotes mean value of variable 



18 


APPENDIX B 

CORE DESIGN RADIAL AND AXIAL TEMPERATURE DISTRIBUTION 

The axial power distribution is assumed to be the same at all radii in 
the core. It has a chopped cosine form as shown in table IV and figure 5, 
with an axial peak to average power ratio of 1.23. The radial peak to aver- 
age power ratio is 1.33. That is, the highest power Zone I fuel pin, the 
centerline pin generates 1.33 times the power of an average pin in the core. 
For the Zone II pin generating the highest power, this ratio is 1.20 and 
for Zone III, it is 1.01. Using these power distributions and the data in 
table I, the curves in figure 5 were obtained showing maximum temperature 
in a cross-section of the fuel as a function of axial location, assuming a 
gap between fuel and clad that is uniform around the pin and a coolant flow 
passage that is uniform around the pin. These indicate that the highest 
fuel temperature in the core, under normal operating conditions, is in the 
centerline fuel pin and slightly downstream of the midpoint of this pin. 

This temperature is normally 1373 K (2006° f)‘. The fuel temperature drops 
off both axially and radially from these values. The maximum fuel tempera- 
ture in Zone II is 1353 K (1982° F) and in Zone III, 1333 K (1950° F) . The 
clad temperatures at the location of maximum fuel temperatures are 1231 K 
(1756° F), 1225 K (1745° F) , and 1215 K (1728° F) for Zones I, II, and III, 
respectively. 

The clad temperature reaches its maximum in each zone at the coolant 
outlet. However, for the calculation of hot spot factors, for both the 
fuel and clad, temperatures at the fuel maximum point were used. It makes 
computations simpler to calculate hot spot factors at the same point for 
the fuel and the clad. It is a valid approach also in the sense that hot 
spot factors for the clad will be conservative at this location, that is, 
greater than calculated at the clad maximum temperature location. 

All the calculations omit any consideration of axial heat transfer. 

We may assume that this has very little effect on results, because a cal- 
culation shows that in the fuel in the region of highest axial tempera- 
ture gradient where there is a heat generation of about 5.9 kW (20 000 Btu / 
hr) to be disposed of, less than 3.5 W (12 Btu/hr) is transferred axially 
out of the fuel. 
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APPENDIX C 

HOT SPOT CALCULATIONS 
Variability of Temperature Rises 

Table V shows details of the calculation of the composite standard de- 
viation for clad temperature and for fuel temperature, using temperatures 
obtained from one dimensional and two dimensional heat transfer calcula- 
tions at the axial position of maximum fuel temperature of the centerline 
fuel pin. Temperature rise from the inlet coolant temperature to the fuel 
hot spot was divided into five increments as shown below. 

1. Coolant temperature rise 

2. Film temperature rise from coolant to clad outer surface 

3. Clad temperature rise from outer surface to inner surface 

4. Temperature rise across the clearance gap between clad and fuel 

5. Temperature rise within the fuel 

Some of the factors contributing uncertainty to the temperatures had 
their effect on all increments, others on only one or two increments. Fac- 
tors affecting the first three items above were used in finding the stand- 
ard deviation for the clad temperature. All 10 factors contributed to the 
standard deviation for the fuel temperature. 

The contribution of each factor (the column headings of table V) to 
the variability of each of the five temperature rise increments was calcu- 
lated as summarized by table V. The calculations used repeated applica- 
tions of equation (3) for each of the five temperature rise increments. 
These five increments exist between , locations that can be defined as fol- 
lows. 


Location 


Subscript, 

j 


Annulus inlet 1 
Local coolant 2 
Clad outer surface 3 
Clad inner surface 4 
Fuel outer surface 5 
Fuel hot spot 6 


The five increments of temperature can also be identified with the 
subscript j as follows: 



Increment 


20 


Subscript , 

j 


Coolant rise 1 
Film rise 2 
Clad rise 3 
Gap rise 4 
Fuel rise 5 


The nominal 
the increment of 


temperature at any location can be 
nominal temperature can be defined 


defined as T m> j and 
as 6 . and thus 



* 

T 

m, j+1 


* 

- T , 
m,J 


(Cl) 


Thus for each of the subtables of table V, the third line of parts (a), 

, V / , , V . . A* A* A« A« 


(b) , (c), (d) , and (e) gives the quantity 
6* respectively 


m,l* m,2’ 


6 „ , 6 . , 
m,3 m,4* 


and 


The next step is to calculate increments of temperature rise <5 m .^ . 
as a function of each independent variable dr factor and the un-’ 

certainty, Ax^, in that variable, from the appropriate functional relation: 


j . . 

m;i,J 


= T . - (x. + Ax, ) - 
m,]+l x i 


m 


,j (x i + 


Ax i ) 


The difference between and <S m \ is the amount by which the 
j rise is changed due to a change in the ith factor by the amount of its 
uncertainty Ax.^. 


The number of standard deviations of x^ asso- 
ciated with Ax-j^ is assumed (as previously discussed) to be n^ (as given 
in table V) and the resulting estimate of the standard deviation, o^, of 
the variable x^ is 


a . 

l 


Ax . /n , 
x x 


(C2) 


Correspondingly, the difference between 6 m . ^ j and <S m • (which had 
been computed from Ax-^) can therefore be reduced £o a value attributable 
to just one standard deviation of x^ by dividing the difference by n^. 
The result is defined as AT m .j_ j: 


AT . . 
m;x,j 


(S , , - S ,)/n, 


(C3) 


Equation (C3) was used to compute the sixth line of each of tables V(a) 
through V(e)o 
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Combination of Variances 


The independent variables were listed as the column headings of 
table V. They are assumed to be statistically independent as one of the 
basic assumptions justifying the use of equation (6). This statistical 
independence is the basis of equation (5) which was specialized to equa- 
tion (6) and is the reason why the increments of equation (6) are added 
together as a sum of their squares. By way of contrast, the temperature 
increments given in the sixth line of each of tables V(a) to (e) for any 
one column are not statistically independent of each other, but are each 
tied functionally to the single random variable given by the column head- 
ing. Because of this nonindependence, their effects are directly additive 
in producing the net temperature effect on the clad or on the fuel. Thus 
the increment of clad temperature variation associated with one standard 
deviation of x^ is the direct sum of the increments given by each of the 
sixth lines (eq. (C3)) of tables V(a) to (c) : 


AT 


cZ ,i 



(C4) 


Results of this summation are given in the first two columns of table III. 


Similarly, the columns headed ATf ^ of table II conform to the same 
basic discussion as do the columns labeled AT C £ ^ of table III. The 
ATf values are merely those summations of the sixth lines of tables V (a) 
through (e) which give the increments of fuel hot spot temperature associ- 
ated with one standard deviation of x^ , as calculated from increments 
given by equation (C3) according to the direct summation: 



(C5) 


The squares (tables II and III) of ATf ^ and AT C £ j, respectively, 
are needed because the next combination of temperature increments will be 
with respect to variables assumed to be statistically independent. Whereas 
equations (C4) and (C5) exhibited summations with respect to j of incre- 
ments that were nonstatistically independent, the combination of those sum- 
mations with respect to i will be over statistically independent variables 
and is done in the manner of equation (5) as specialized to equation (6) . 

In particular, equation (6) was specialized to the variance of the clad 
inner temperature as equation (C6) and to the variance of the fuel hot spot 
temperature as equation (C7) 


10 



i=l 


2 

°T y cSL 


(C6) 
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10 


T , f 


r£ 

: i=l 


< 1T f,i y 


(C7) 


The results of the computations with equations (C6) and (C7) on the 
entries in tables III and II, respectively, are given in the bottom lines 
of the same tables. The standard .deviations (which are the square roots 


of the variances a 
tables III and II. 


T,c£ 


and 


aj f) are given in the footnotes to 


Individual Hot Spot Factors 


A definition of an individual hot spot factor was given as equa- 
tion (4). To perform an approximate worst case or "deterministic" analysis, 
the numerator of equation (4) is replaced by the i^ term of equation (3) 
with the Ax£ in that term replaced by 3a^. Under these assumptions, the 
numerator of equation (4) is given by 


3 T 

9x7” U°i) - 341 
1 7 


(C8) 


where ATf was given in table II. Correspondingly, "individual hot spot 
factors" were calculated as given in table VI from temperature increments 
&Tf^ given in table II, according to the equation 


f . 

i 


3AT f,i 

(T ■ - T ) 

fm coolant inlet'ref 


= 1 + 


(C9) 
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TABLE .1. - FUEL PIN HEAT TRANSFER DATA 
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TABLE II. - FUEL HOT SPOT DATA 


Contributing variable 

AT 


(AT f 

,i> 2 

K 

°F 

K 2 

V 

1. Overall power adjustment 

8.3 

— 

68.9 

223.2 

2. Local power 

4.2 


17.6 

57.2 

3. Local fuel thickness 

0.6 

1 

0.36 

1.17 

4. Local thermal conductivity 





in fuel 

6.5 

11.7 

42.3 

136.9 

5. Orifice mismatch 

1.3 

2.3 

1.69 

5.5 

6. Flow annulus width 

5.6 

10.1 

31.4 

101.6 

7. Heat transfer coefficient 

0.45 

0.81 

0.20 

0.66 

8. Fuel-clad gap width 

28.2 

50.8 

795.2 

2577 

9. Fuel shift in clad 

2.0 

3.6 

4.0 

12.9 

10. Clad shift in flow annulus 

1.3 

2.3 

1.69 

5.5 

2 

a m 



963.3 

3121.2 

T>f 






a T>f = 31.0 K (55.8 F) 
Hot spot factor =,1 + 


3a, 






coolant inlet 


) 


= 1.46 


ref 


fm 


coolant inlet 


ref 1370 K (2006° F) 


1166.5 K (1640 F) 
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TABLE III. - CLAD HOT SPOT DATA 


Contributing parameter 

AT cJt,i 

'«T cl ,l> 2 

— 

K 

— 

°F 

K 2 

°F 2 

1. Overall power adjustment 

2.7 

4.9 

7.3 

23.6 

2. Local power 

1.4 

2.5 

2.0 

6.4 

3. Local fuel thickness 

— 




4. Local thermal conductivity 





in fuel 

— 




5. Orifice mismatch 

1.3 

2.3 

1.7 

5.5 

6. Flow annulus width 

5.6 

10.1 

31.4 

101.6 

7. Heat transfer coefficient 

0.45 

0.81 , 

0.20 

0,66 

8 . Fuel-clad gap width 

— 




9. Fuel shift in clad 

6.4 

11.5 

41.0 

132.7 

10. Clad shift in flow annulus 

8.9 

16.0 

79.2 

256.6 




162.7 

527.5 


a TjCjl 12.8 K (23.0° F) 


Hot spot factor 

= 14- 

(T - T ) 

elm coolant inlet re £ 

T clm 

ref 1231 K (1760° F) 

> 

T 

coolant inlet 

1166.5 K (1640° F) 


1.59 
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TABLE IV. - AXIAL POWER DISTRIBUTION 


Relative distance 
from inlet end 

Power relative to 
average pin power 

0.0 (inlet) 

0.605 

O.i 

.815 

.2 

.990 

.3 

1.3.2 

.4 

1.20 

.5 

1.23 

. 6 

1.20 

.7 

1.12 

.8 

.990 

.9 

.815 

1.0 (outlet) 

.605 







TABLE V. - CALCULATION OF NOMINAL AND PERTURBED TEMPERATURE RISES, 


>» 

I— I M 4J 
« 0» -H 
O 3 CO 

.o o c 

hJ a qj 


to a 4J 3 
a ft o 
o a> 3 

J ^ T) C 
U C 
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Calculated uncertainty 




TABLE V. - Concluded. CALCULATION OF NOMINAL AND PERTURBED TEMPERATURE RISES 



h ft.in ft 

a) 6 

3 <U* 0 0 

Cm +J •o *o 


































TABLE VI. - "DETERMINISTIC" (APPROX. WORST CASE) TEMPERATURE 

AND HOT SPOT FACTOR 


Contributing variable 

Individual hot spot factor, 
f . 

l 

1. Overall power adjustment 

1.122 

2. Local power 

1.062 

3. Local fuel thickness 

1.009 

4. Local thermal conductivity in fuel 

1.096 

5. Orifice mismatch 

1.019 

6. Flow annulus width 

1.083 

7. Heat transfer coefficient 

1.007 

8. Fuel-clad gap width 

1.416 

9. Fuel shift in clad 

1.030 

10. Clad shift in flow annulus 

1.019 


10 

TT f. = total hot spot factor = 2.176 
i=l 1 

Maximum fuel temperature = 1610° F 
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FIGURE 3 - NODE. ARRANGEMENT FOR 
SYMMETRIC CASES 


Figure 4- node arrangement for 

asymmetric cases 
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